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The current responsivity of a semiconductor superlattice THz-photon detector is calculated using 
an equivalent circuit model which takes into account the finite matching efficiency between a detector 
antenna and the superlattice in the presence of parasitic losses. Calculations performed for currently 
available superlattice diodes show that both the magnitudes and the roll-off frequencies of the 
responsivity are strongly influenced by an excitation of hybrid plasma-Bloch oscillations which are 
found to be eigenmodes of the system in the THz-frequency band. The expected room temperature 
values of the responsivity ( 2-3 A/W in the 1-3 THz-frequency band) range up to several percents of 
the quantum efficiency e/hu of an ideal superconductor tunnel junction detector. Properly designed 
semiconductor superlattice detectors may thus demonstrate better room temperature THz-photon 
responsivity than conventional Schottky junction devices. 
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I. INTRODUCTION 



Theoretical investigations of high-frequency properties of semiconductor, superlattices have been carried out for 
almost thirty years starting from the trail-blazing works by Esaki and TsuEJfl. They proposed to prepare superlattices 
made of different semiconductor materials in order to realise an artificial periodical system that would allow to observe 
Bloch oscillations due to Bragg reflection of electrons from the boundaries of the Brillouin zone. Making use of the 
semiclassical Boltzmann equation they found that the dc current- voltage characteristics of the superlattice should show 
, a negative differential conductance in sufficiently strong electric fields, namely, when the Bloch frequency Qb = eEd/h 
(where e is the electron's charge, E is the dc electric field in the superlattice n d is the superlattice period, h is the 
Plank's constant) becomes comparable to the effective scattering frequency mII. They proposed to make use of the 
ON ■ negative differential conductance for development of an oscillator in the THz-frequency bandEl and also suggested that 
the superlattice might posses promising, properties as an artificial material for non- linear mixing of photonscl. 

Later, several papers have proposecHcj to probe Bloch oscillationS|-jn a superlattice in a strong dc electric field 
by applying an additional strong ac field, r-^clf- induced transparency^, dc current suppression by a strong ac field 
causing dynamical localisatiop,.of carriersfjtij, absolute negative conductance^, and current harmonics generation with 
oscillating power-dependenceaB have been predicted to occur in the presence of resonant interaction between Bloch 
Q | oscillating electrons and a strong THz-frequency radiation. 

A profound analogy between the-dynamics of electron wave-packets in periodic systems and the dynamics of super- 
conductkig Josephson junctionst3'e3 was put forward to explain the main features of the superlattice THz- frequency 
responseEJ. It has been suggested, in particular, that the resonant steps on the dc current- voltage characteristics of 
the irradiated superlatticesE-J should occur in analogy with 'Shapiro steps' normally observable in irradiated super- 
conducting junctionstfl. Several recent papers show that the curious behaviour of electrons in the superlatticjes may 
give rise to a rich variety of new non-linear phenomena occurring in the THz- field irradiated superlatticescH~Ea ._. ._. 

In a number of papers Esaki- Tsu .negative differential conductance was found to be experimentally realisableE§E£l 
and has been carefully examineoOcJ in the context of development of new millimetre-wave band (0.03-0.3 THz) 
oscillators. On the other hand, the recent observations (jf^a strong dc current suppression indicating dynamical 
localisation of electronsE3, absolute negative conductanceEiH23, and 'Shapiro steps' on the dc current-voltage curve of 
the THz-field irradiated superlatticesES open the prospects for applications of the superlattices as novel solid state 
detectors operating in 1-10 THz-fcequency band which the Bloch frequency in the superlattices normally belongs too. 

It has recently been estimateoE3 that the room temperature current responsivity of a superlattice detector ideally 
coupled to the THz-photons can nearly reach the quantum efficiency e/fiu (where u) is the incident radiation frequency) 
in the limit of high frequencies u > v. This value of the responsivity is being normally considered as a quantum 
limit for detectors based on superconducting tunnel junctions operating at low temperaturesEEI. For high frequencies 
the mechanism of the THz-photons detection in superlattices was describecE3 as a bulk superlattice effect caused by 
dynamical localisation of electrons. 

In the present paper we develop a self-consistent theory of the superlattice current responsivity. We apply the Boltz- 
mann equation approach for describing the electron motion in the superlattice minibando and assume an equivalent 
circuit for the superlattice coupled to a broad-band antenna (see Fig. [l]), which is similar to the equivalent circuit used 
in resonant-tunnellingEa and Schottky diodescH simulations. The suggested equivalent circuit of the device allows one 



1 



to treat microscopically the high-frequency response of the miniband electrons and, simultaneously, take into account 
a finite matching efficiency between the detector antenna and the superlattice in the presence of parasitic losses. Ops 
analytic results lead to the identification of an important physical concept: the excitation of hybrid plasma-BlochH 
oscillations in the region of positive differential conductance of the superlattice. The numerical computations, per- 
formed for room temperature behaviour of currently available superlattice diodes, show that both the magnitudes and 
the roll-off frequencies of the responsivity are strongly influenced by this effect. The excitation of the plasma-Bloch 
oscillations gives rise to a resonant-like dependence of the responsivity on the incident radiation frequency, improving 
essentially the coupling of the superlattice to the detector antenna. We will also show that peak current densities 
in the device and its geometrical dimensions should be properly optimised in order to get maximum responsivity for 
each frequency of the incident photons. Finally, we will present numerical estimates of the responsivity for the 1-4 
THz frequency band and compare its value with the quantum efficiency e/fiu> of an ideal detector. 



II. DESCRIPTION OF THE MODEL 



For the description of the ac electr|On-,transport in a superlattice we use a quasi- classical wave packet treatment of 
the electron motion in a superlatticetrEI. The energy spectrum of electrons in a miniband is taken in a tight-binding 



approximation: 
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where A is the superlattice miniband width, d is the superlattice period, p z is the quasi-momentum of an electron along 
the superlattice axis (perpendicular to the layers) , p x ,Py are the quasi-momentum components along the superlattice 
layers, and m is the effective mass of electrons along the superlattice layers. 

The quasi-classical velocity v z (p z ) of an electron moving along the superlattice axis and the time-derivative of the 
quasi-momentum are given by the expressions 

v z (jp z ) = —r — = Wo sin — - (2) 



dp z \ h 

Pz = eE z (t) (3) 

where e is the electron charge, vo = Ad/2h~ is the maximum velocity of electrons along the superlattice axis, and E z (t) 
is the time-periodic electric field directed along the superlattice axis. 

The electric current density j z (t) is calculated from the distribution function, 

m-efvMfM-frjfc, (4) 
which satisfies the time-dependent Boltzmann equation 

Below we use the relaxation-time approximation for the collision integra. 

'df\ /(p,*)-/o(p) 



at 



(6) 
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where r is the constant relaxation time for electron's scattering, and fo(p) is the equilibrium distribution function. 
The relaxation-time approximation is, of course, an oversimplification of the numerous scattering processes taking 



place in a real superlattice. Nevertheless, several paperscfl have demonstrated that the phenomenon of dynamical 
localisation can be described very well within this approximation when compared to the corresponding results obtained 
from a full-scale Monte-Carlo simulation. This circumstance lends support to the simplified model for the collision 
integral. The great advantage of Eqs.(||) and (||) is that they allow an analytical calculation of the time-dependent 
current, to be used in the equivalent circuit analysis, and thus we can study in detail various parametric dependencies 
of the calculated quantities. 

We would also like to emphasise her&,tiiat the wave packet description of electron motion in a superlattice is valid 
if the following inequalities are fulfilled™!! 
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eEd <C A, Aq (7) 
?iw<A,A G (8) 
?u/<A,A G (9) 

were A G is the width of the superlattice minigap, v = 1/r is the electron's collision frequency, and u> is the charac- 
teristic frequency of the external ac field. We also assume that the electronic mean free path £pp = vqt = Ad/2Tiv 
is smaller than the superlattice length L in order to neglect the influence of the boundaries on the superlattice high- 
frequency properties. 

III. THEORETICAL FORMALISM 
A. Path integral 

The exact .solution of Eqs. (||) and (||) for arbitrary time-dependent electric field can be presented in the form of a 
path integralii 

/(P,*)= / i/dtiexp [-i/(t -ti)]/ (p*- f eE z {t 2 )dt 2 ) (10) 

Using Eqs.(Q) and (|Io| ) we find the time-dependent current I(t) describing ac transport in a superlattice, with electron 
performing ballistic motion in a mini-band according to the acceleration theorem and suffering scatteringB 



/t r r-t 

vdt x exp [-u(t- h)] sin — / V(t 2 )dt 2 
-oo |_iV/l J tl 



(11) 



where V(t) = LE z (t) is the voltage across the superlattice perpendicular to the layers, L = Nd is the superlattice 
length, N is the number of periods in the superlattice sample, I p — Sj p , S = ira 2 is the area of the superlattice, a is 
the superlattice mesa radius, and 

v f ( Pzd\ . , , 2dp 

is the characteristic current density. The integration over p z in Eq. ( |l2| ) must be carried out over the Brillouin zone 
—irh/d < p z < irk/d. n 

In particular, for a dc voltage Vsl applied to the superlattice (EgjlT)) results in the Esaki-Tsru current- voltage 



curve 



t sl (v \ 0T (Vsl/Vp) 
1dc\ v sl) = 21 P — ——2 (13) 
1 + {Vsl/VpY 

At peak voltage Vsl — Vp = Nfiv/e or, equivalently, at peak electric field E$l = Ep = hv/ed, the current in the 
superlattice reaches its maximum (peak) value Ip , so that jp can be defined as the peak current density. 

The peak current density jp and the scattering frequency v can be considered as the main parameters of the 
employed model. They can readily be estimated from experimentally measured or numerically simulated values of Ip 
and Vp. For both degenerate and non-degenerate electron gas one getsBB 

3p = eny (14) 

if A > kT : €F, where kT is the equilibrium thermal excitation energy, ep = h 2 (3n 2 n) 2 / 3 /(2m e g) is the Fermi 
energy of degenerate electrons, n = J / (p)2dp / (2Trh) 3 is the charge carrier density, m c g = rri'J z Z m}l' i is the density 
of states effective mass near the miniband bottom, and m zz = 2h 2 /Ad 2 is the effective mass of electrons along 
the superlattice axis. In the particular case of the Boltzmann equilibrium distribution function Eq.(|l^) yieldscl 
jp = (envo/2)[Ii(A/2kT)/I (A/2kT)], where Jo,i are the modified Bessel functions. 
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B. Monochromatic excitation 



We now suppose that in addition to the dc voltage Vsl, an alternating sinusoidal voltage with a complex amplitude 
Vaj is applied to the superlattice: 



V(t) = V SL + \ K exp(»wt) + V* oq>(-iwt)] 



(15) 



Generally, Vsl , V w can be found from an analysis of the equivalent circuit given in Fig. ^. We write the ac voltage 
amplitude as V u = |K,|e 1 ^; both \yj\ and ip can be obtained self-consistently taking account of reflection of the THz 
photons from the superlattice and their absorption in the series resistor Rg. 
Making use Eq. (|ll|) we obtairfl: 



I(t) = 2I P / vdti exp 
Jo 



(—vti) sin 



eV S L 



Nfi 



ti + $(*,ti) 



where 



e 1 

$(Mi) = x „ {iV u exp{iu}t)[exp(-iu}ti) - 1] + c.c.} 

Alfa 2 



(16) 



(17) 



According to Eq.(|l6|), electrons in a superlattice miniband perform damped Bloch oscillations with the frequency 
Qb = eVsL/Nh — eEsLd/U, and the phase <I>(i, ti) modulated by the external ac voltage. 

In the limit of V u — > Eq.(^6|) reduces to the dc current-voltage curve given by Eq. (|l3|) . On the other hand, when 
the frequency of the ac voltage is small, wt < 1, we get 



eti 1 

$0, h) = x 2 [ V " exp(i^i) + c.c] 



(18) 



which corresponds to a slow modulation of the Bloch frequency by external voltage. In this case the current in the 
superlattice follows instantaneously the time- dependent ac voltage according to the dc current-voltage curve. 

Equation ( |l6| ) contains, as special cases the following results: (i) a harmonic voltage ^(t) (Vsl — 0) leads to 
dynamical localisation, and current harmonics generation with oscillating power dependencetl; (ii) a dc current- voltage 
characteristics of the irradiated superlattice. f&c (Vsl, VJ) — (ui/2tt) J I(t)dt shows resonance features ('Shapiro steps') 
leading to absolute negative conductanceEHj'EZl; (iii) and to generation of dc voltages (per one superlattice period) that 
are multiples of hcu/S3. 

C. Method of perturbations 

Let us assume that the external ac voltage Vj is so small that perturbation theory holds, while the dc voltage Vsl 
applied to the superlattice keeps its finite value. Expanding Eq. (16) around Vsl at V u — * in a Taylor series, we 
obtain the time-dependent electric current in the form: 



Kt) = 1 S dc(Vsl) + \ [G S A U", Vsl)V^ 1 + c.c] + A/g£(w, V SL ) 



where 



G ! &(u > Vsl) = G F 1 (cj,Vsl) 



is the superlattice ac conductance^, Go = 2Ip/Vp is the superlattice conductance at u) — > 0, Vsl — * 0, 

1 + iLOT- (Vsl/Vp) 2 



Fi(uj,V S l) 



[1 + (Vsl/Vp) 2 } [(1 + iurf + (V SL /V P y 



(19) 



(20) 



(21) 



is a dimensionless function describing the dependence of the superlattice conductance both on frequency and applied 
dc voltage (field), and 



MffrfaVsL) = -M 2 F 2 (uj,Vsl) 



(22) 
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where 

F , v , _ I s d l c (Vsl + Nhw/e) - 2I S D L C (V SL ) + I s d l c (Vsl - Nhu/e) 
2{ ' SL> {NUio/eY 

_ Alp (Vsl/Vp) [3 + (lot) 2 - (Vsl/Vp) 2 ] 

Vl [1 + {V S l/VpY\[1 + {Vsl/Vp + lot) 2 ][1 + (V SL /V P - lot) 2 } 

is the change in the dc current in a superlattice caused by THz-photonsS. 
At u> — > Eq.(20) yields the dc differential conductance of the superlattice 



(23) 



1 - {Vsl/Vp) 



= Go v r ' 2 (24) 

[1 + {Vsl/Vp) 2 } 2 

while at Vsx — » one gets the well known Drude formula for the ac conductivity of the electron gas <j(u>) 

a(w) = (25) 

1 + ILOT 

where the small-field dc conductivity of the superlattice ao can be given in terms of the dc conductance gq — GqL/S 
or, equivalently, in terms of peak current density and the peak electric field 

a = 2g (26) 

According to Eqs.(|l5|) and ( |l9| ) the ac power Pj^ absorbed in the superlattice is 

= \^G S ac^, V SL )K\ 2 (27) 

whereptlie real part of the superlattice conductance can be presented in analogy with Eq. (23) in a finite difference 
formcreS 

ReG^( W) V S l) = [ iS dc(Vsl + NfiLo/e) - I S D L C (V SL ~ NfiLo/e)] 



2NHlo 

2Ip 1 + (lot) 2 - (Vsl/Vp) 2 

V P [1 + {Vsl/Vp + lot) 2 }[1 + {V SL /V P - lot)' 



(28) 



We note here that our calculations assume a uniform dc/ac electric fields inside the superlattice. This assumption is 
justified for the sub-threshold dc electric fields Vsl < Vp, andior for ac field frequencies larger than the transit-time 
frequency of space-charge waves u> > vq/L in the superlatticeSa. 



D. Stimulated emission and positive photo-current 

Equation (^8) for the superlattice conductance demonstrates a possibility of the resonant interaction of the THz- 
photons with electrons performing damped Bloch oscillations in the superlattice. For lot 3> 1 the condition for the 
resonance has a form lo — ±(Vsl/Vp)t which can be equivalently expressed as to = ±Q_b- At dc electric fields 
Esl > Epy/l + {lot) 2 the real part of the superlattice conductance starts to be negative which corresponds to 
negative absorption (stimulated emission) of the THz-photons. In the limit lot 1 photons with frequency lo > Qb 
are absorbed and the ones having frequency lo < 0# are emitted. 

Equation (p2|) describes the dc current change in the superlattice under the influence of the THz-photons. The 



standard result of the classical rectification theorylia 



AI s D L c (co -> 0,V SL ) = \\V^\ 2 I s D L c {Vsl) (29) 



where 



Ip {Vsl/Vp) [3 - (Vsl/Vp) 2 } 



'dc ( v sl) - -4=j 3 

Vp l + (V S l V P y\ 



(30) 
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is the second derivative of the dc current- voltage curve of the superlattice, is obtained from Eq. (|22) when the electron 
system relaxes during one period of the ac field. 

At dc bias fields E$l < Ep^/Z + {lot) 2 the current change in the irradiated superlattice is .negative which cor- 
responds to the onset of suppression of the current due to dynamical localisation of carriersOTLj. In the opposite 
case Esl > Ep^/3 + {lot) 2 a positive current change (positive photocurrent) occurs. It is im po rtant to note that 
the condition of occurrence of the positive photo-current implies, according to Eqs.(p2|) and (p7|), the existence of 
stimulated emission of photons from the superlattice, P^ s < 0. 

The increase of the dc-current in superlattices caused by stimulated emission of THz-photons has recently been 
observed in experimentsEl These experimental results lend support to the validity of the theoretical formalism 
employed in the present paper for the analysis of the superlattice THz-field response. 



IV. CURRENT RESPONSIVITY 

A. Equivalent circuit calculation 

We define the current responsivityB of the superlattice detector as the current change AI induced in the external 
dc circuit per incoming ac signal power P,. 

Ri(u,V SL ) = ^j- ( 31 ) 

This definition takes into account both the parasitic losses in the detector and the finite efficiency for impedance- 
matching of the incoming signal into the superlattice diode. In the further analysis of the current responsivity we use 
both the frequency / = tojlix and the angular frequency to notations. 

According to Eqs.^2] and ^7], in the small-signal approximation both the dc current change A/j^ and the power P^ s 
absorbed in the superlattice are proportional to the square modulus of the complex voltage |V^,| 2 . This circumstance 
permits us to calculate self-consistently |Kj| 2 for given values of the incoming power making use a linear ac equivalent 
circuit analysis and, then, find the current responsivity Ri{u), Vsl)- 

The results of the calculation of the superlattice current responsivity Ri(ui, Vsl) are presented in the following form: 



Rf^VsiM^Vsi 
Rs(dI^ L c {VsL)/dV S L) 



rmv sl ) = /\^t::^\ (32) 



where 



,(0), xr s e {Vsl/Vp)[3 + {lot) 2 - {V SL /Vp) 



Nhv [1 + (V SL /V P ) 2 }[1 + {lot) 2 - {Vsl/Vp) 

is .the superlattice current responsivity under conditions of a perfect matching and neglecting parasitic losses {Rs — > 
0)El ,where R s is the series resistance. 

The factor A{lo,Vsl) in Eq. ( |32| ) describes the effect of the electrodynamical mismatch between the antenna and 
the superlattice and the signal absorption in the series resistance 

Z a -{Z^{lo,V sl ) + R s ) 2 



A{lo,V S l) = 



1 - 



Z A + {Z s Mu,V SL )+Rs) 



ReZ s A ^V SL ) 
ReZSL(u,V SL ) + Rs ' 1 j 



The first factor in Eq. (f54|) describes the reflection of the THz-photons due to mismatch of the antenna impedance 
Za and the total impedance of the device Vsl) + Rs, with the second one being responsible for sharing of 

the absorbed power between the active part of the device described by the impedance Z^{lo,Vsl) and the series 
resistance R$- 

The superlattice impedance is defined as 

Z^{lo, Vsl) = 1/ [G s a l c {lo, V S l) + iwC] (35) 

where G A c{lo, Vsl) is the superlattice conductance, C — eoS/AnL is the capacitance of the superlattice, and eo is the 
average dielectric lattice constant. 

Finally, the last factor in Eq. ([52j) describes the redistribution of the external bias voltage Vnc between the dc 
differential resistance of the superlattice {^^{Vsl) / cLVsl) -1 and the series resistance -Rft^with the dc voltage drop 
on the superlattice Vsl being determined by the solution of the well-known load equatiorH 

Vdc = Vsl + Idc{V S l)Rs ■ (36) 



G 



B. Classical rectification at high frequencies 



Suppose now that the following set of inequalities is satisfied 

LOT < 1 (37) 

Rs « Ro (38) 
uj » 1/y/RoRsC , (39) 

which implies that the frequency of the incident radiation u> is small in comparison with the scattering frequency of 
electrons is, series resistance Rs is small in comparison with the small-field dc resistance of the superlattice Ro = 1/Go 
and the appropriate RC time of the device is high enough. In this case Eq. ( |32|) gives 

V SL ) = 2I^\V SL) x {1+ Z z A a/Rs)2 x jj^y 2 (40) 

Equation (^) was obtained by Sollner et.alH following Torrey's and Whitmer's approach!]] developed for analyses 
of classical crystal rectifiers. They also took account of the mismatch between the antenna impedance Za and the 
device impedance. It was applied to analyse detection of photons by resonant tunnelling diodes in the THz- frequency 
range. This expression demonstrates that the responsivity of the classical rectifier at high frequencies is proportional 
to the second derivative of the dc current-voltage curve, and hence strongly decreases with increasing of frequency, 
Ri{uJ,V SL ) oc uj- 2 . 

The factor 1/RsC in Eq. jib] ) defines the roll-off frequency for the current responsivity and, consequently, capac- 
itance of the device should be minimised in order to increase the responsivity in the high frequency range. The fast 
decrease of the responsivity with increasing of frequency can be attributed to imperfect matching of the device to the 
antenna impedance when incoming photons either are reflected from the device or are absorbed in the series resistor. 

In our case Eq. ( f40| ) shows that the current responsivity of the superlattice can reach its maximum (negative) 
value at bias voltage Vsl — OAVp, i.e. in the region of the positive conductance. Further, the responsivity of the 
superlattice does not depend on the length of the device because Vp oc L and C oc 1/L, and, finally, the responsivity 
is proportional to the current density jp. Consequently, in order to increase the responsivity in the high-frequency 
range one should employ highly conductive superlattice samples for which inequalities (37)- (39) may be not satisfied. 
In this case the interaction of the superlattice with the incoming radiation can be essentially changed by excitation 
of the eigenmodes in the superlattice device. 



C. Superlattice dielectric function. Hybridisation of Bloch and plasma oscillations 

Let us analyse the condition of optimised matching of the superlattice to the incident radiation going back to the 
general Eq. (J3^). Assuming the limit of negligible series resistance R,s — » this condition can be obtained from 
solution of the equation 

Z^(uj,V S l) = Z A (41) 

for the complex frequency a; (Vsl) ■ This solution determines the resonant line position and the line width at which 
the absorption in the superlattice tends to its maximum value. 

Using Eqs. ( pof) and (|3^) one can transform Eq. ( [4l| ) to the following form: 

e(u J ,E SL ) = ^— (42) 



where 



e{u,E SL ) = e + ^F 1 {io,E SL ) (43) 



is the dielectric function of the superlattice, with the dc field E$l being applied to the deviceEa, Fi(u>, Esl) is defined 
by Eq. ©. 

In the high-frequency limit cq/CZalo — > the solution of Eq. ( |4l| ) coincides with the solution of the equation 

e(w, E SL ) = (44) 
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describing eigenfrequencies w± of the hybrid plasma-Bloch oscillations in a superlatticd^l 

1/2 



SL) 



1 + {Esl/Ep) 



— ) (Esl/Ep) 2 

Up J 



(45) 



where cup is the plasma frequency of electrons in a superlattice. The plasma frequency cup can be given in terms of 
the small-field dc conductivity ctq or, equivalently, in terms of the peak current density jp 



UJp 



( 4it<jQv\ 



Snjped 
e h 



1/2 



(46) 



Equation ( |4q ) reduces in the particular case of wide-miniband superlattices (A 3> kT,ep) to the standard formula 
u) P = (4Tre I n/e m zz ) 1 /' 2 . 



In the limiting case of small applied dc electric fields Esl/Ep — ■+ one finds from Eq. (45) the plasma frequency 



u>± — > ±Wp, while in the opposite case Esl/Ep — > oo, the Bloch frequency — > ±f^ = ±eEsLd/Ti is recovered. 
The scattering frequency ^ in Eq. (0) is responsible for the line width of the plasma-Bloch resonance. 

We have calculated the hybrid plasma-Bloch oscillation frequency fjj = u^;/2ir, using Eqs. @ and @, for the 
typical values of the superlattice parametersE3~cJ eo — 13, d ~ 50A, Ep ~ 10 kV/cm, = ^/27r = 1.2THz for 
different values of the current densities j p (Fig. ||). For small values of the current densities jp ~ 10 kA/cm 2 the 
frequency of the hybrid oscillation increases with applied voltage in all range of the parameter Vsl/Vp ■ On the other 
hand, for higher values of the current densities jp ~ (50 — 1000) kA/cm 2 the hybrid oscillation's frequency starts to 
decrease with increasing bias voltage in the sub-threshold voltage range Vsl < Vp. Then, at super-threshold voltages 
Vsl > Vp, luh starts to increase again tending to the Bloch frequency. It is important to note that at high values 
of the dc current densities jp the hybrid plasma-Bloch oscillations become well defined eigenmodes of the system 
ifn > fv)- Therefore, an essential improvement of the matching efficiency between antenna and the superlattice can 
be expected in the high- frequency range due to a resonant excitation of this eigenmode in the device. 



D. High-frequency limit 

In the high-frequency case, when the signal frequency / — lu/2it is larger than the scattering frequency /„ = v/2w, 
description of the superlattice response based on the classical rectification theory is no longer valid. Let us employ 
here Eq. (|3^) in order to analyse the high-frequency limit of the responsivity in the ideal case of the perfect matching 
and neglecting the series resistance. At u 3> i.e. in the frequency band were ac field absorption and negative 

photo-current are predominated in the superlattice response, Eq. ([33|) yields the frequency- independent value for the 
current responsivity 

d(°)c \r \ e Vsl/Vp 

R] '{UJ -> OC,VsL = -T77- -I T 7T7 /T/ x 2 47 

Nhv 1 + (Vsl/VpY 

At positive bias VsL/Xif-t> the responsivity is negative (reduction of the dc current occurs due to dynamical 
localisation of carriersatLa). The dc voltage dependence of the responsivity reproduces the dc Esaki-Tsu current- 
voltage curveS. The responsivity reaches its maximum value i^max = — e/(2-/V7jz/) = — l/(2Vp) at Vsl — Vp. For 
high values of the applied voltages Vsl ^> Vp we get r\° '(u — > oo, Vsl — * oo) = — 1/Vsi,- 

Let us compare the high-frequency limit of the responsivity of the superlattice with the quantum efficiency rJ? max = 
e/hw which is believed to be a fundamental restriction for the responsivity of superconductor tunnel junctionals. This 
quantum efficiency, (or quantum limit) corresponds to the tunnelling of one electron across the junction for each signal 
photon absorbedliS, with a positive sign of the responsivity. 

In our case the mechanism of the photon detection is different (see Fig. ||) . Electrons move against the applied dc 
electric force due to absorption of photons. At Vsl — Vp the responsivity is negative, indicating that one electron is 
subtracted from the dc current flowing through the superlattice when the energy 2eVp is absorbed from the external 
ac field. One half of this energy is needed for the electron to overcome the potential barrier which is formed by the 
dc force, with another half being delivered to the lattice due to energy dissipation. If the applied dc voltage is strong 
enough, i.e. Vsl 3> Vp, dissipation plays no essential role in the superlattice responsivity. In this case the energy 
eVsL should be absorbed from the ac field in order to subtract one electron from the dc current simply due to the 
energy conservation law. 
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In order to demonstrate typical frequency scales involved in the problem we plot in Fig. Q the frequency dependence 
of the responsivity -R\°\u>, Vsl) for a superlattice with d = 50A, N = 40, L = 0.2/xm, E P = lkV/cm, A = 20meV, 
fv = v /2ir = 0.12THz at sub-threshold voltage Vsl = 0.99Vp as a function of the frequency / = ui/2w. This example 
roughly corresponds to the superlattice samples experimentally investigated in Ref. |33J. The responsivity decreases in 
the region / < f v and, then, tends to the constant value determined by Eq. (]47|). At frequency Jql — 2Nv/2ti the 
value of the responsivity equals the quantum efficiency e/hf. 

However, the semiclassical approach employed in the present paper is restricted by the inequality / < /a — A//i. 
This inequality simply requires that there must be an allowed transition between different Wannier-Stark states in 
the superlattice miniband due to photon absorption as illustrated in Fig. || It is important to note that /a//ql = 
£fp/L < 1 (where £pp — v t = Ad/2hi> is the electron's mean free path) in our calculations and, therefore, the current 
responsivity of the whole sample is always less than the quantum efficiency e/hf. In a wide-miniband superlattice 
with A — > 2Nhv (or, equivalently, in a short superlattice sample with L — > tpp ) the responsivity is approaching the 
quantum efficiency in the high-frequency limit. 

As was mentioned above, an increase of the length L of a detector may have an advantage in the reduction of 
the parasitic capacitance C oc 1/L, and, hence, in the increasing of the roll-off frequency X/RsC of the device. It is 
interesting to note that in case of N superconductor tunnel junctions connected in series the resulting responsivity 
is expected to be e/iVftwjnplying that one photon should have been absorbed in each junction in order to add one 
electron to the dc circumLa. In the case of a superlattice samples consisting of N unit cells (periods) the resulting 
responsivity can be presented as —e/2NTiv = (—e/NTuS) x (lo/2v). For uj > 2v the responsivity of the unit superlattice 
cell starts to be larger than the quantum efficiency e/huj. This conclusion can be readily understood if one realises that 
due to absorption of one photon in the miniband electrons can 'jump' over several superlattice periods as illustrated 
in Fig. | 



V. RESULTS AND DISCUSSION 



In this section we shall investigate the current responsivity of the superlattice making use of the Eq. ([52|) which takes 
into account both the finite matching efficiency between the antenna and the superlattice and the parasitic losses in 
the superlattice diode. We present our results using the dimensional value of Ri(w, Vsl) (given in units of Ampere / 
Watt) and also introducing the normalised responsivity 

R tN (u, V SL ) = Ri(u, V SL )/(e/Tuj) (48) 

which permits us to compare directly the responsivity R^uj^Vsl) of the superlattice with the quantum efficiency 
e/Tiuj . In oik calculations we use the typical geometrical dimensions of the superlattice samples investigated in the 
experimeijitsP'Tcj and also assume that the bow-tie antenna impedance — 50f2 does not depend on THz-photon 
frequency^. 



A. Excitation of the plasma-Bloch oscillations 



For demonstration of the frequency dependence of the superlattice current responsivity in the THz-frequency band 
we will focus on the GaAs/Gao.sAlo.sAs superlattices specially designed to operate as millimetre wave oscillators at 
room temperature. In Ref. ^ wide-miniband superlattice samples with d = 50A, A ~ 113meV, n ~ 10 17 cm -3 , were 
investigated experimentally. They demonstrated a well-pronounced Esaki- Tsu negative differential conductance for 
Esl > Ep ~ 4kV/cm with the high peak current of the order of jp ~ 130kA/cm 2 . The measured value of the peak 
current is in a good agreement with the estimate jp ~ (80 — 160) kA/cm 2 for n ~ (1 — 2) x 10 17 , T = 300K based 
on Eq. (|l^), if one assumes an equilibrium Boltzmann distribution for the charge carriers. From the peak electric 
field and current we find the scattering and plasma frequencies j v ~ 0.5THz, fp = 2THz, respectively, assuming 
eo = 13 for the average dielectric lattice constant. The maximum frequency for the semiclassical approach to be valid 
for these samples is /a — 27 THz. Figure shows the frequency dependence of the normalised current responsivity 
calculated for three values of the peak current density in the superlattice, i.e. jp = 13, 130, and 300 kA/cm 2 and for 
three values of the peak electric field, Ep = 4, 9, and 13 kV/cm. We also use the typical values for the superlaitice 
length L = 0.5 (im (superlattice consists of 100 periods), and assume a = 2/im for the superlattice mesa radiusCJ l3. 
We choose Rs = 10fi for the series resistance of the device ia the THz- frequency band, i.e. the same value as for 
resonant tunnelling diodes having the same radius of mesasE3. The calculations are performed in the region of the 
positive differential conductance for dc bias voltage close to the peak voltage (Vsl = 0.95Vp) 
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For Ep = 4kV/cm (/„ ~ 0.5 THz) Fig. || demonstrates well-pronounced resonant behaviour of the normalised 
responsivity as a function of frequency. The resonance frequency and the maximum value of the responsivity rise if 
the peak current density increases. For j p = 300 kA/cm 2 the normalised responsivity reaches its maximum value 
—RiN — 0.02 (—Ri ~ 2 A/W) at frequency / ~ 2.5 THz . For higher values of the peak electric fields Ep = 9 
kV/cm (/„ ~ 1.08 THz) and Ep = 13 kV/cm (f v ~ 1.57 THz) the resonance line-widths are broadened due to 
implicit increase of the scattering frequencies. In particular, for Ep — 13 kV/cm, jp = 300 kA/cm 2 the normalised 
responsivity has an almost constant value —RiN — 0.006 (—Ri ~ 0.6 A/W) up to / ~ 2.5 THz and, then, rapidly 
decreases. 

The frequency behaviour of the normalised responsivity originates from excitation of the plasma- Bloch oscillations 
in the superlattice. We indicate in Fig. |H] the positions of the hybrid frequencies fjj = \w± \/2ir with arrows. For small 
peak electric fields (low values of the scattering frequencies) the hybrid frequency corresponds to the maximum of the 
normalised responsivity For higher values of the peak field (higher values of the scattering frequencies) it corresponds 
to the roll-off frequency at which the responsivity starts to decline. 

Frequency dependence of the responsivity for different applied dc fields , Vsl = 0.3, 0.6 and 0.95 Vp is illustrated 
in Fig. ^ for jp = 130kA/cm 2 . The same geometrical dimensions and circuit parameters of the superlattice device 
have been chosen for calculations as previously. The normalised responsivity decreases with decreasing of the bias 
field tending to zero at Vsl — * 0. On the other hand, the position of maximum responsivity shifts to lower values 
with increasing of the bias field in full agreement with Eq. (|45|) . 

B. Optimised superlattice length 

The enhancement of the normalised responsivity requires an optimum matching efficiency of the superlattice to 
the broad-band antenna and minimisation of the parasitic losses in the series resistor. These requirements impose an 
optimum length of the superlattice for each chosen frequency of the incoming THz-photons and series resistance. 

We show in Fig. the dependence of the normalised responsivity on the number of the superlattice periods for 
/ = 2.5THz.We used for calculation a = 2/im, jp = 130 kA/cm 2 , Vsl = 0.95Vp, and three values of the series 
resistance Rs = 10, 30 and 50 ft. For all three values of the series resistance the responsivity displays a well 
pronounced maximum for the optimum number of the superlattice periods N — N max . The value of N max increases 
with increasing of the series resistance ( N max ~ 40 for Rs — lOfi, N max ~ 60 for Rs = 30fi and for iV max ~ 90 for 
Rs = 50fi). This result can be readily understood by recalling that a larger volume of the superlattice minimises 
parasitic losses for higher values of the series resistance because of reduction of the sample's capacitance. 

For incoming photon higher frequencies the parasitic losses in the superlattice device play even more important role 
leading to a further increase of -/V max . Figure |S| shows the dependence of the normalised responsivity on the number 
of the superlattice periods for / = 3.9THz and the same superlattice parameters as in Fig. [?]. In this case (iV max ~ 70 
for R s = 10O, iV max ~ 110 for R s = 30fi and N mgcK ~ 180 for R s = 50Q ). We can conclude, therefore, that the bulk 
mechanism of the superlattice high-frequency response provides important benefits for operation of the superlattice 
detectors in the THz-frequency band. 

C. Responsivity dependence on bias voltage 

Experimental studies of interaction of high-frequency fields with the miniband electrons in semiconductor superlat- 
tices having a relatively small (jp ~ 10 kA/cm 2 ) current densities were performedocj for frequencies / = 90 GHz, 
/ = 450 GHz and / = 3.5 THz. A field-induced reduction of the current through the superlattice was reported, 
with a different bias dependence below and above a characteristic frequency (ITHz). Distinct bias dependence of the 
responsivity was attributed to quasistatic (ujt <C 1) and dynamic (lot 3> 1) interaction of the miniband electrons with 
ac fields at frequencies below ancL-above ITHz, respectively. The intraminiband relaxation time (r ~ 10~ 13 s) was 
estimated from experimental dataei 

In Fig. H we plot the dc bias dependence of the responsivity at / = 450GHz and / = 3.9THz for the superlattice 
parameters corresponding to experimentsca (jp ~ 10 kA/cm 2 , Ep ~ 13 kV/cm, a ~ 2.5/um, L ^i-Q.5//m, Rs — 20S1) 
together with the dc Esaki-Tsu current-voltage curve. In a full agreement with the observationsc£l, at / = 450 GHz 
responsivity reaches its maximum value at Vsl — 0.5Vp (roughly corresponding to the position of the maximum of 
the second derivative of the dc current-voltage curve), while for 3.9 THz-radiation the position of the maximum is 
shifted to the peak voltage Vp . This can be readily understood if one takes into account that for ujt > 1 the model of 
the classical rectification is no longer valid. In this case the current change under THz-photon irradiation is described 
by Eq. (p2[) giving the second derivative of the current in the finite difference form taking account of finite photon 
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energies . At high frequencies uit — > oo this equation yields AI^(u>, Vsl) oc I^c^Vsl), i.e. the bias-field dependence 
of the responsivity should reproduce the dc current-voltage curve showing the maximum value at Vsl — Vp. 

It is important to note that with increasing peak current densities this behaviour qualitatively changes. Figure [h] 
demonstrates the bias-field dependence of the responsivity for the same frequencies and superlattice parameters as 
Fig. ^ and for jp ~ 130kA/cm 2 . First, we emphasise that the responsivity is considerably higher than in the previous 
case both for 450 and 3.9THz radiation. Second, for 450 GHz-radiation the responsivity reaches its maximum value 
at Vsl — Vp and not at Vsl — 0.5Vp as previously. For the high peak current densities the superlattice impedance 
variation due to applied dc voltage essentially changes the coupling efficiency between antenna and superlattice. This 
leads to a qualitatively different behaviour of the responsivity for low and high peak current densities: the latter can 
manifest in experiments a more efficient coupling of radiation into the superlattice. 



D. Optimised peak current density 

For a given superlattice geometry the current responsivity of the superlattice can also be enhanced by choosing an 



optimum value of the peak current densities. This circumstance is illustrated in Fig. 11 where we plot the responsivity 
as a function of the peak current density for Vsl — 0.99Vp, Ep ~ 4 kV/cm, a ~ 2.5/^m, L ~ 0.5/xm, Rs — lOfi and 
for three frequencies of radiation, i.e. / = 1, 2 and 3 THz. For different frequencies the responsivity reaches its peaks 
at different values of j™ ax (/) which increases with increasing of the radiation frequency. 

This behaviour can also be explained by excitation of the plasma-Bloch oscillations in the superlattice by THz- 
photons if one takes into account the resonant condition / = /p(jp ax ), where the hybrid frequency fn = \u±/2ir\ 
is given by Eq. (||) (see also Fig. |). Calculated values jp^lTHz) ~ 60kA/cm 2 , j P lax (2THz) ~ 220kA/cm 2 and 
jp ax (3THz) ~ 540kA/cm 2 are shown in Fig. [ll] by arrows. They are in a good agreement with the positions of the 
peaks found from calculations based on the general Eq. (j32|). 

We can conclude, therefore, that high current density superlattices should be used in order to achieve the large 
responsivity values (2-3 A/W) in the 1-3 THz -frequency band. The high current densities can be obtained by 
choosing pffijide-miniband and/or highly-doped samples. For example, according to Eq. (Jl4|) in superlattices with A ~ 
130 meVc3 the peak current density jp ~ 1000 kA/cm 2 can be reached for n ~ 3 x 10 17 cm" 3 . In this case the 
equilibrium thermal excitation energy (kT ~ 26 meV) and the Fermi energy of degenerate electrons (ep ~ 27 meV) 



are considerably smaller than the miniband width. Hence, Eq. (14) can be employed for the peak current density 
estimates both for room and low-temperature conditions. 

Finally we would like to note that the measurements of the responsivity of resonant-tunnelling (double- barrier) 
heterostructure diodes reported in Ref. [|3^| were carried out at frequencies ras high as 2.5 THz. The reported value 
of the responsivity was of the order of several /LtA/W. As has been mentioned! 3 ^, this value is smaller by over an order 
of magnitude than the responsivity of THz-frcquency Schottky diodes in this frequency band. We believe that an 
optimised superlattice detector as discussed above may have some advantages over these devices due to rather high 
expected current responsivity in the THz-frequency band. 



VI. CONCLUSIONS 



In conclusion, we have calculated the current responsivity of a semiconductor superlattice detector in the THz- 
frequency band based on equivalent circuit modelling. Using a path integral solution of the Boltzmann equation 
within the relaxation-time approximation for the collision integral we have obtained an analytical expression for the 
responsivity, taking account of (i) frequency dependent superlattice response to a THz field caused by dynamical 
localisation of electrons, (ii) a finite matching efficiency between a detector antenna and the superlattice, and (iii) the 
presence of parasitic losses in the device caused by a series (contact) resistance. 

We find that the responsivity of the superlattice ideally coupled to the incident radiation tends to a finite value 
with increasing radiation frequency. This value is simply determined by the energy conservation law governing the 
THz-photons absorption in the presence of scattering processes. 

Excitation of plasma-Bloch oscillations in the superlattice (which arc found to represent eigcnmodcs of the system in 
the THz-frequency band) can essentially enhance both the magnitude and the roll-off frequencies of the responsivity 
due to resonant coupling of radiation into the superlattice. The excitation of the plasma-Bloch oscillations can 
manifest itself as a resonance-like dependence of the normalised responsivity on the THz- photon frequency and (or) 
as a specific dependence of the responsivity on bias fields for superlattices showing high peak current densities. 

Changes in peak current density and in superlattice length can effect drastically the coupling efficiency and parasitic 
losses in the superlattice device. These parameters are found to play an important role in the THz- field detector 
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performance and need to be optimised for each value of the radiation frequency and series resistance. For higher 
frequencies the optimum length of the superlattice detector tends to increase demonstrating a benefit of a bulk-type 
mechanism of the THz-photon detection. 

In currently available superlattices possessing the optimised length and the peak current the responsivity is expected 
to be as high as ~ (2-3) A/W in the 1-3 THz-frequcncy band. These values of the responsivity range up 
to several percents of the quantum efficiency e/hcu of an ideal superconductor tunnel junction for this frequency 
range. The analysis of the current responsivity performed in the present paper does not necessarily assume cooling 
of the superlattice: the estimated values of the current responsivity of superlattices can be expected at even room 
temperature. 
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FIG. 1. (a) THz-radiation coupled to a iV-period semiconductor superlattice by a co-planar broad band bow- tie antenna, 
Pi and P r are the incident and reflected powers respectively, (b) Equivalent circuit for a THz-photon detector with a dc 
voltage bias source: _B-miniband electrons capable to perform Bloch oscillations, C-superlattice capacitance, i?g-parasitic 
series resistance, Za -bow-tie antenna impedance, Vd c -dc bias voltage. 
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FIG. 2. The calculated hybrid plasma-Bloch oscillation frequency fn is plotted as a function of the normalised superlattice 
voltage drop Vsl/Vp for different values of the peak current densities jp = 10, 50, 100, 300, 500, and 1000 kA/cm 2 . Typical 
values of the superlattice parameters ( d — 50A, _Ep=10 kV/cm, eo =13) were used for the calculations. 
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FIG. 3. Real space energy diagram illustrating THz- photon (/ ^> v/2-k) detection in the superlattice: DC electric field Esl 
is applied to the iV-period semiconductor superlattice with the miniband width A . Under the action of the dc field electrons 
perform Bloch oscillations with the spatial amplitude A/eEsL ■ At critical dc electric voltage (field) Vsl — Vp — NTiv/e 
(eEshd = hv) electrons move against the dc electric force due to absorption of photons climbing up the Wannier- Stark ladder. 
The energy 2eVp should be absorbed from external ac field in order to subtract one electron from the external circuit. One 
half of this energy is needed for the electron to overcome the potential barrier which is formed by the dc force, with the other 
half being delivered to the lattice due to energy dissipation. A quasi-classical description of the process is valid if / <C A/h 
when allowed transitions between different Wannier-Stark state exist. 
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FIG. 4. The calculated current responsivity —R ( f \f) of the ideally coupled superlattice with d — 50 A, N = 40, Ep — 
1 kV/cm, A= 20 meV at sub-threshold voltage Vsl = 0.99Vp (solid curve), and the quantum efficiency e/hf (dashed curve) 
are plotted as a function of the frequency / = lu/2h. At frequency /ql = 2Nv/2-n the value of the responsivity equals the 
quantum efficiency. The calculation of the responsivity is valid up to the frequency /a = A/ft which should be smaller than 
fq L = 2Nu/2-k within the framework of the employed model. In a wide-miniband superlattice A — > 2N x hp the responsivity 
can approach the quantum efficiency in the high-frequency limit. 
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FIG. 5. The frequency dependence of the normalised current responsivity \Rin\ = \Ri/(e/hui) \ of the superlattice THz-photon 
detector (a = 2/im, L = 0.5/im, Rs = 1057, Vsl = 0.95Vp) is calculated for three values of the peak current density (jp = 
13, 30, and 300 kA / cm 2 ) and for three values of the peak electric field (Ep — 4, 9, and 13 kV/cm ). The relevant positions 
of the hybrid plasma-Bloch frequencies fn are indicated for each curve by arrows showing characteristic resonance (high peak 
current densities) and roll-off (low peak current densities) behaviour. 
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FIG. 6. The frequency dependence of the normalised current responsivity \RiN\ = \Ri/(e/huj)\ of the superlattice THz-photon 
detector (a — 2/im, L = 0.5/im, _Rg = 10f2, _Ep =4kV/cm, Jp = 130 kA/cm 2 ) is calculated for three values of the applied dc 
voltage Vsl — 0.3, 0.6, and 0.95 Vp. The hybrid plasma-Bloch frequencies fn (indicated for each dc voltage by arrows) are 
found to be in a good agreement with positions of maximum values of the responsivity. 




FIG. 7. The dependence of the normalised current responsivity \IUn\ = \Ri/{ e /^ luJ ) \ of the superlattice THz-photon detector 
(a = 2/im, Ep = 4kV/cm, jp = 130 kA/cm 2 ) is calculated at / = 2.5 THz for three values of the series resistance (Rs=W, 30, 
and 50 fi) as a function of number of the superlattice periods N. 
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FIG. 8. The dependence of the normalised current responsivity \RiN | = \Ri / (e/hui) | of the superlattice THz-photon detector 
(a = 2/im, Ep = 4kV/cm, jp = 130 kA/cm ) is calculated at / = 3.9 THz for three values of the series resistance (Rs=W, 30, 
and 50 fl) as a function of number of the superlattice periods N. 
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FIG. 9. (a) Normalised Esaki-Tsu I — V characteristics of a superlattice THz-photon detector experimentally investigated 
in Ref. 41 (a — 2/im, Ep — 13kV/cm, jp = 10 kA/cm 2 , L — 0.5/im, Rs = 20f2). Current responsivity is calculated for this 
detector for (b) 450-GHz and (c) 3.9-THz radiation as a function of the normalised dc voltage Vsl/Vp. 
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FIG. 10. (a) Normalised Esaki-Tsu I—V characteristics of a high current density superlattice THz-photon detector (a = 2/im, 
Ep = 4 kV/cm, jp — 130 kA/cm 2 , L = 0.5/xm, Rs = 2057). Current responsivity is calculated for this detector for (b) 450-GHz 
and (c) 3.9-THz radiation as a function of the normalised dc voltage Vsl/Vp. 
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FIG. 11. The current responsivity of the superlattice THz-photon detector (a = 2.5/im, Ep = 4.5kV/cm, L — 0.5/im, 
Rs = 20£1) is calculated as a function of the peak current density for three frequencies of radiation: / = 1, 2and 3 THz. The 
relevant positions of the peak current densities j™ ax calculated from criterion of excitation of the plasma-Bloch oscillations 
/ = /n(jp ax ) are shown for each frequency / respectively. 
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